Generation of large-amplitude coherent-state superposition via ancilla-assisted 

photon-subtraction 
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We propose and demonstrate a novel method to generate a large-amplitude coherent- state su- 
perposition (CSS) via ancilla-assisted photon-subtraction. The ancillary mode induces quantum 
interference of indistinguishable processes in an extended space, widening the controllability of 
quantum superposition at the conditional output. We demonstrate this by a simple time-separated 
two-photon subtraction from continuous wave squeezed light. We observe the largest CSS of trav- 
elling light ever reported without correcting any imperfections, which will enable various quantum 
information applications with CSS states. 

PACS numbers: 42.50.Dv, 42. 50. Ex, 03.67.-a, 03.65.Wj 
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Superposition of macroscopically distinct coherent 
states, often called coherent-state superposition (CSS), 
is regarded as a realization of Schrodinger's famous cat 
paradox. Typical CSS states are defined as \C±) = 
-J==(\a) ± |— a)), where | ± a) are coherent states with 

amplitudes ±a. The |C+) (|C_)) is called an even- (odd- 
) CSS state and \a\ 2 is often regarded as its "size" since 
it reflects the distance of two superposed coherent states. 
Such states have been realized in a few physical systems 
P, Q. Among them, large size CSS states of travelling 
light are important for many quantum information tasks 
such as linear-optics quantum computation 0, Q , quan- 
tum teleportation [5[, and quantum metrology [q|. In 
practice, one can generate approximate \C±) by subtract- 
ing photons from a squeezed vacuum p}. Along this line, 
single-photon subtraction has been demonstrated, gener- 
ating odd cat-like states [1, 0, Qij • But their sizes ware 
still small, typically \a\ 2 ~ 1.0. 

Although subtracting more than two photons leads to 
larger states, experiments get more challenging due to 
rapid decrease of the success probability. Recently an 
alternative way utilizing a photon-number state and ho- 
modyne detection was proposed and a generation of the 
state y/2/3 |2) — \/l/3 |0) from a pulsed two-photon state 
was demonstrated [H|. It is close to a superposition of 
two squeezed states where they are displaced from the 
origin in opposite directions by the amount ~ 1.2 in the 
phase space. If 3.5 dB squeezing is further applied, the 
state would be |C+) with \a\ 2 ^ 2.6. However, it remains 
a challenge to apply squeezing onto such a squeezed CSS 
state. 

In this letter, we propose and demonstrate a novel way 
to enhance the size of CSS states without resorting to fur- 
ther subtraction of more than two photons or squeezing 



operations. Instead we introduce an ancillary mode to 
assist in suppressing the weights of smaller number pho- 
tons in a CSS state, and hence to enhance its size. In 
the experiment, we implement it in the time domain, by 
time-separated two-photon subtraction from continuous 
wave (cw) squeezed light. 

Photon subtraction is done by tapping a small fraction 
of the squeezed vacuum for photon counting, and by se- 
lecting the transmitted state conditioned on the detection 
of photons. A two-photon subtraction from a single-mode 
squeezed vacuum (without ancillae) is described as 



a 2 S(e )\0) =(3 S(e ) ( 
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where a is an annihilation operator describing one-photon 
subtraction, S is a squeezing operator, eo represents the 
degree of squeezing [12| and (3q = eo/(l — eft). This is a 
squeezed state of the superposition V2eo |2) + |0), which 
well approximates an even CSS state with the size up to 
\a\ 2 ~ 1.0 [7]. However, it is still possible to further in- 
crease the size if one could optimize the ratio between |0) 
and 1 2) in the superposition, independently from en [131] . 

Suppose two independent photon-subtraction events 
instantaneously occur at time t\ and respectively, on 
a cw squeezed vacuum generated from an optical para- 
metric oscillator (OPO) as depicted in Fig. Ufa). Since 
the cw squeezed vacuum has a finite bandwidth £ , the 
effect of each photon subtraction spreads over the tem- 
poral wavepacket ij^t — t\.2) m = V^ e_ ^°^ _£l ' 2 ' m the 
transmitted beam 
tion A = |^2 — h 
ets are highly overlapped and two temporal waveforms 
defined as *±(t) = [ip(t - h) ± ip(t - t 2 )]/y / 2(l ± I A ) 
with Ia = (1 + CoA)e~^° A appear as orthonormal mode 
functions to describe the process [13, Oil . The symmet- 
ric mode acts as the main mode to be measured 



14 [TBI . fl6j . When the time separa- 
satisfies (qA < 1, these two pack- 
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FIG. 2: (color online) Theoretical Wigner functions of the 
CSS state generated by the time-separated two-photon sub- 
traction and their photon number distributions, e = 0.3. (a) 
CoA = and (b) CoA = 1.4 correspond to approximate |C+) 
of |a| 2 = 1.2 and 2.6, respectively. 



FIG. 1: (color online) (a) The time-separated two-photon sub- 
traction and the temporal mode functions \£±(t). OPO: op- 
tical parametric oscillator, APD: avalanche photo diode, (b) 
Schematic of the ancilla- assisted photon subtraction, (c) A 
dependence of the squeezing e + (solid line) and e_ (dashed 
line), e = 0.3. 



while the asymmetric mode Sfr-(t) serves as an ancilla 
and the essential mechanism can be well described by 
these two modes. Precisely speaking, subtracted photons 
have finite correlations with the other temporal modes 
in the transmitted beam, but the effects of those modes 
can simply be taken as a small optical loss [l3[. Then 
the temporal two-mode model can be translated into a 
simplified spatial two-mode model with modes *f?± and 
a 50/50 beam splitter (BS) as illustrated in Fig. [TJb). 
The coincidence of clicks at each photon detector always 
means that two photons must come from either or 
and not that single photon from each mode. This 
is due to the bunching nature of photons. These two 
subtraction processes are indistinguishable, producing a 
superposition, 

(a 2 + -a 2 _)S + (e + )S.(e.)\0) + \0)_, (2) 

where the subscripts ± indicate modes \£±. The e± rep- 
resent effective degrees of squeezing in each mode which 
are determined by a fully multimode theory and depend 
on A, Co, and the OPO pumping parameter e |l3[. Fig- 
ure [T](c) shows typical behaviours of e± as functions of 
A. Since e_ is relatively small at Co A < 1, using Eq. (pQ), 
Eq. ([2|) can be rewritten as follows by neglecting terms 
proportional to e 2 _ or higher orders, 

Eq. © « 

(3+S + (e + ) (e+a 1 , 2 + (l - ^Jj |0> + S_(e_) |0)_ 
^A/-|$) + £_(e_)|0>_, (3) 



where (3+ = e+/(l — e+), a normalized state |$) and 
a normalization factor M are introduced. Compared 
with Eq. ([1]), it can be seen that the ancilla mode pro- 
vides additional controllability for tuning the superpo- 
sition in the main mode through e_, which is simply 
controlled by changing the time separation A as shown 
in Fig. [U(c). This allows us to coherently modify the 
photon number distribution of the state in the main 
mode and in particular suppress the weight of small pho- 
ton number components to increase the size of the CSS 
state. Note that, in a rigorous cw model, the state in 
the main mode is slightly degraded to a statistical mix- 
ture p+ = (l- Co)|*)(*| + C 5+|0)(0|5| even without 
practical imperfections due to a weak entanglement be- 
tween the main and ancilla mode. Figure [2] compares the 
theoretical Wigner functions of p + based on the full cw 
model [13[ for ( A = and 1.4, i.e. without and with the 
ancilla assistance. These correspond to the CSS states of 
H 2 = 1.2 (CoA = 0) and 2.6 (CoA = 1.4) with the fi- 
delities of 0.928 and 0.932, respectively, where the latter 
A is chosen to maximize the size of CSS while preserv- 
ing a high fidelity (> 0.9). For larger A, Co becomes 
non- negligible as discussed later. 

A Schematic of our experimental setup is shown in 
Fig. [3j It is built on our previous single-photon sub- 
traction setup [10[ with qualitative updates to cope with 
long data acquisition for hours. A cw Ti: Sapphire laser 
of 860nm is used as a light source. A squeezed beam is 
generated from an OPO which employs a periodically- 
poled KTiOP0 4 (PPKTP) as nonlinear crystal. The 
OPO is pumped by a frequency doubled beam of 20mW 
from a second harmonic generation (SHG) cavity with 
a KNb03 crystal. This corresponds to e = 0.3. The 
OPO bandwidth is Co/2tt -4.5 MHz (HWHM). A small 
fraction (10%) of the beam is tapped by a tapping BS 
(TBS) and going through two filtering Fabry-Perot cav- 
ities of 2mm- long and of 0.9mm- long, respectively. The 
unwanted photons in non-degenerate modes extending 
over a wide range of frequencies are well filtered out by 
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FIG. 3: (color online) Experimental setup of the time- 
separated two-photon subtraction. SHG : second harmonic 
generation cavity, OPO : optical parametric oscillator, AOM 
: acousto-optic modulator, FC : fiber coupler, TBS : tapping 
beam splitter (R =10%) FBS : 50/50 fiber beam splitter, 
and APDs : Si avalanche photo diodes (Perkin-Elmer SPCM- 
AQR-14). 



these cavities and only the component in the degenerate 
mode around 860nm is guided into two avalanche photo 
diodes (APDs) 

The beam transmitted through the TBS is continu- 
ously measured by a homodyne detector. All cavities are 
actively locked on the resonance for 860nm by a weak 
coherent beam and electronic feedback. The coherent 
beam is periodically chopped by an acousto-optic modu- 
lator (AOM), which defines distinct time bins by presence 
and absence of the coherent beam, and thus enables the 
alternate sequence of locking and photon counting with- 
out stray light. Also via phase locking the coherent beam 
to parametric gain of the OPO, we obtain phase infor- 
mation of the measured state from interference between 
the coherent beam and the local oscillator (LO) beam at 
the homodyne detector. 

Clicks of APDs announce that a conditional state ap- 
pears in the transmitted beam. The time separation A 
is adjusted by an electronic delay line. The pumping of 
e = 0.3 at the OPO outputs a 3.0dB squeezed vacuum in 
temporal mode ^+(t) and the simultaneous-click rate of 
the two APDs is about 6 counts per second for A = 0ns 
and less than 1 count per second for A ~ 100ns. The 
conditional state in the main mode is measured by a ho- 
modyne detector to reconstruct its Wigner function. For 
every pair of clicks of the APDs with desired A, a contin- 
uous photo-current of the homodyne detector is stored. A 
quadrature amplitude in the main mode is extracted from 
the continuous data by applying the numerical mode fil- 
ter with the shape of ^+(t) which is equivalent to an 
optical LO with the waveform of ^+(t) Hol|. 

Experimentally reconstructed Wigner functions and 
photon number distributions of p + with A = and 32ns 
are given in the upper row of Fig. [U(a). The recon- 



struction was performed by the conventional maximum 
likelihood estimation method [13] from 22,000 points of 
quadrature amplitudes in the raw experimental data. No 
experimental imperfections are corrected in the recon- 
struction procedure, whereas most of the previous works 
corrected the losses at the state evaluation step [H, . 
Figure[D(a) clearly shows that with finite A, the two posi- 
tive peaks corresponding to the components | ±a) become 
apart, i.e. the size-increase of a. 

Ideally, the two-photon subtracted states of A =0 and 
32ns (Co A = and 0.9, respectively) would approximate 
|C+) of \a\ 2 =1.2 and 2.0, respectively, with high fi- 
delities (> 90%). In the experiment, the fidelities of 
76% and 60%, for \a\ 2 =1.2 and 2.0, respectively, are 
achieved. A particular feature of the macroscopic quan- 
tum superposition is the presence of negative dips in the 
Wigner function due to the interference between two dis- 
tinct wavepackets. For A =32 ns, the negative dips are 
clearly observed at the sides of the origin. This is, to 
our knowledge, the largest CSS state of the travelling 
light ever observed, showing the negative Wigner func- 
tion without any compensation of experimental imper- 
fections. Such quality is achieved by use of a highly pur e 
squeezed vacuum generated from our low loss OPO |lo|. 
The sources of imperfection are separately investigated. 
The OPO has 0.4% intra-cavity loss. Associated with 
9.5% transmittance of the coupling mirror, it gives 0.95 
of the escape efficiency. The propagation efficiency from 
the OPO to the homodyne detector is 0.95. The quan- 
tum efficiency of the homodyne detector and its elec- 
tronic noise give effective detection efficiency 0.98. The 
spatial overlap between the squeezed vacuum and the LO 
gives another efficiency of 0.96. In total, the overall ef- 
ficiency is estimated around 0.85. Theoretical Wigner 
functions including these imperfections are shown in the 
lower row of Fig. S^a). The experimental results are al- 
most in agreement with the theoretical models [18] while 
the negativity is slightly degraded by some uncleared im- 
perfections. We also note that the maximum CSS state 
is obtained around A = 48ns (see Fig. [2](b)) while the 
experimentally observed negative dips are limited up to 
32ns since the low count rates at large A also degrades 
the negativity. 

The photon number distribution in Fig. SJa) also 
reveals a mechanism of the ancilla-assisted photon- 
subtraction. Suppression of the small photon-number 
components by increasing A is clearly observed. Our 
result should also be contrasted to the squeezed state 
superposition ^/2/3\2} — ^/T/3\0) generated in pjj. In 
our case, the even- number distribution extends at least 
up to n = 4, and is clearly deviated from that of the 
squeezed vacuum, which should be another signature of 
the macroscopic nature of the superposition. 

FigureHJb) plots the average photon numbers (n±) of 
modes \£± as a function of A, that agree with the theory 
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FIG. 4: (color online) (a) Wigner functions of the state p+ in mode ^+ for A = Ons (CoA = 0, left) and A = 32ns (CoA = 0.90, 
right). Upper row: experimental Wigner functions (W(x,p)), their contours, and photon number distributions (p(n)). Lower 
row: model calculation of the Wigner functions, (b) Average photon numbers in modes \£±. Upper (filled) plots: experimental 
(n+) with A = 0, 32, 48, 64, 96 ns, from left to right. Lower (unfilled) plots: experimental (n_) with A = 32, 48, 64, 96 ns, 
from left to right. Solid and dashed lines are the theoretical predictions including 15% optical loss. 



well indicated by lines [l8[. The increase of (n+) in the 
range (oA < 2 supports an experimental evidence of the 
size- increase. In this range most of the photons accumu- 
late in the main mode with fewer photons in the ancillary 
mode. As A increases beyond CoA > 2, however, (n_) 
starts to increase, and photons in modes and \£_ are 
getting entangled. It induces degradation of the purity 
of the reduced state in the main mode. For (oA ^> 1, 
the state in the main mode is degraded to a completely 
mixed state p + \S + (\2){2\ + 1 0> <0 1 ) 

Finally, we point out that the above drawback of the 
purity degradation could be overcome if one would take 
advantage of an ancillary mode which is not entangled 
to the main mode. This is done by preparing a coherent 
state | a) in the ancillary mode. The improved scheme is 
represented in terms of generic two modes A and B as 

(4 - a%)S A \0) A \a) B = {a\ - a 2 )S A \0) A \a) B . (4) 

Since the coherent state is an eigenstate of a, one can 
make a superposition of two- and zero-photon subtrac- 
tions in mode A without entangling two modes, i.e. with- 
out degradation. Here the coherent amplitude simply 
plays the role of A in the time-separated two-photon 
subtraction. In practice, the experiment will get more 
involved since one has to lock the relative phase of these 
two modes. Also a related scheme was proposed recently 

0. 



In conclusion, we have proposed a concept of the 
ancilla-assisted photon-subtraction and experimentally 
demonstrated it by the time-separated two-photon sub- 
traction from cw squeezed light. Due to the quantum 
interference assisted by the ancillary mode, we have suc- 
cessfully observed and characterized the largest CSS state 
of travelling light so far showing the negative dips of the 
Wigner function without any corrections of measurement 
imperfections. Our experimental scheme involves the 
trade-off between the size-enlargement and the degrada- 
tion of the state purity. We pointed out that the trade-off 
can be circumvented by extending it to the one with a 
coherent-state ancilla. 
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